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Abstract 

Background: Limit dextrinase inliibitor (LDI) inliibits starcli degradation in barley grains during nnalting because it 
binds with limit dextrinase (LD). There is a wide genetic variation in LDI synthesis and inactivation during barley 
grain development and germination. However, the genetic control of LDI activity remains little understood. 

Results: In this study, association analysis was performed on 162 Tibetan wild accessions by using LDI activity, 835 
Diversity Arrays Technology (DArT) markers and single nucleotide polymorphisms (SNPs) of the gene HvLDI encoding 
LDI. Two DArT markers, bpb-8347, bpb-0068, and 31 SNPs of HvLDI were significantly associated with LDI activity, 
explaining 10.0%, 6.6% and 13.4% of phenotypic variation, respectively. Bpb-8347 is located on chromosome 6H, near 
the locus of HvLDI, and bpb-0068 is located on 3H. 

Conclusions: The current results confirmed the locus of the gene controlling LDI activity and identified a new DArT 
markers associated with LDI activity. The SNPs associated with LDI activity may provide a new insight into the genetic 
variation of LDI activity in barley grains. 

Keywords: Limit dextrinase inhibitor (LDI), Genome-wide association study (GWAS), Single nucleotide polymorphism 
(SNP), Tibetan wild barley 



Background 

Malting is the first step in brewing or beer production 
of malt barley, and the properties of malting barley are 
generally evaluated on the basis of malting quality. High 
malting quality is characterized by high malt extract 
(ME) and diastatic power (DP), as well as high wort 
fermentation by yeast, i.e. complete degradation of 
starch and polysaccharide in malt [1]. The major chemical 
compound in barley grains is starch, which consists of 
30% amylose and 70% amylopectin [2]. 

Starch is degraded by coordinative action of a- amylase, 
|3-amylase, limit dextrinase (LD) and a-glucosidase in the 
endosperm of the germinating grains [3]. LD is the only 
enzyme of cleaving a- 1-6 linkages in branched dextrins 
molecules [4]. Without the action of LD, the branched 
dextrins could not be fermented by yeast [5]. Hence, LD 
plays an important role in the brewing and distilling 
industries for producing fermentable sugar. 
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LD in barley grains exists in three different forms: 
insoluble bound, soluble inactive and active free, and 
only the active free form participates in starch mobilization 
[6]. However, the free LD accounts for a small portion of 
total LD. The low activity of available LD in grains is due 
to its combination with endogenous inhibitors, limit 
dextrinase inhibitor (LDI) [7]. LD and LDI interact with 
each other in a 1:1 molar ratio complex [8]. The LD 
bound to LDI is thought to be a limiting factor for 
complete degradation of starch [9]. 

LDI is synthesized in the developing grains, and degraded 
gradually during the malting process [7]. Consequently, LD 
activity in malt is not only dependent on the amount of the 
enzyme itself in mature grains, and also is closely related 
to the amount and type of LDI. HvLDI was identified to 
encode the LDI protein, and located at chromosome 6H 
[10,11]. However, genetic variation and controlling of LDI 
synthesis and degradation in barley grains and malt is not 
totally understood. In view of a wider variation of LD 
activity in Tibetan wild barley than in cultivated barley 
(unpublished data), it is imperative to determine the 
genetic variation of LDI content in Tibetan wild barley. 
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which has been recently proved to be a progenitor of 
cultivated barley and is rich in genetic diversity [12,13]. 

Association analysis has been widely used to identify 
genes or loci of important traits, such as abiotic stress 
tolerance [14,12,15], agronomic traits [16,17] and also 
barley quality [18,19]. However, there has been no report 
on the discovery of novel loci or elite alleles attributed 
to LDI activity using association analysis. 

The objectives of this study were to (1) determine the 
genetic variation of LDI activity in Tibetan wild barley; 
(2) identify the DArT markers and SNPs of HvLDI associ- 
ated with LDI activity in Tibetan wild barley. 

Methods 

Plant materials 

A total of 162 Tibetan wild barley accessions, kindly 
provided by professor Sun of Huazhong Agricultural 
University, China, were used for LDI activity and asso- 
ciation analysis. All accessions were planted in early 
November 2011 and 2012 in Zijingang campus, Zhejiang 
University (Hangzhou, China). Every genotype was culti- 
vated in accordance with local agronomic practices with 
three replications. The 162 barley accessions were planted 
in a block with each accession consisting of three lines 
(2 m length and line distance is 0.25 m). 

Assay of LDI activity 

Grain samples were micro-malted in a Phoenix System 
Micro-malting Apparatus (Adelaide, Australia) in the order 
of steeping, germination and drying. LD was partially puri- 
fied from barley malt according to Kristensen et al. [20]. 
The fractions were applied to gel filtration chromatography 
after the ion exchange chromatography step, and partially 
purified LD was collected and used for the measurement of 
LDI activity. 

LDI activity was determined according to MacGregor 
et al. [21] with some modification. One ml of 0.1 M sodium 
acetate (pH 5.5) containing 10 mM 1,10-phenanthroline 
was added to 0.1 g barley powder, incubated at 4°C for 
30 min. The extract was heated at 70°C for 40 min, 
centrifuged and the supernatant was collected. The 
protein content of the extracts was measured using 
Bradford assay kits (Sangon Biotech). Twenty micrograms 
protein of LDI extract were mixed with 10 mU partially 
purified LD and the volume was made up to 0.5 ml in 
0.1 M maleic acid containing 0.02% Na azide (pH 5.5). 
The remaining LD was determined using the Limit- 
Dextrinase assay kit (Megazyme). The LDI activity was 
calculated as the reduced LD activity. 

Population structure and kinship analysis 

Totally, 771 DArT markers, with minor allele frequency 
(MAF) higher than 0.03 (Additional file 1: Table SI), were 
used for population structure analysis using STRUCTURE 



software (v2.3.3) [22], in which the number of clusters (k) 
was set from 1 to 12 and ten iterations were performed in 
an admixture model with 10,000 burning period and 
100,000 MCMC (Markov Chain Monte Carlo). The DArT 
markers used were derived from Diversity Arrays Tech- 
nology Pty Ltd, Australia and distributed over the whole 
genome [23-25]. The most probable number of clusters 
(k) was estimated according to the value of Ak. When Ak 
had the highest value, the value of k was the number of 
clusters [26]. Kinship was estimated using SPAGeDi soft- 
ware [27]. Genetic distance and neighbor-joining tree were 
developed with NTSYSpc (version 2.10e). 

Variance component and heritability estimation 

Variance analysis was conducted using SPSS software. 
The used model was: y = mu + GENOTYPE + ERROR 
(fixed in low case, random in capitals). Heritability was 
estimated according to = Vg/(Vg + Ve). 

Genome-wide association study (GWAS) of LDI activity 

Genome-wide association study between LDI activity 
(mean value of 2011 and 2012) and DArT markers was 
conducted using TASSEL software (v3.0), where Q, K 
and Q + K methods were applied [28] . The structure 
matrix was included as covariate to correct population 
structure in Q model: j = Qv + + e, where X is the 
DArT marker matrix, Q is the structure matrix, ^ and v 
are coefficient vectors for DArT marker vector and 
population structure vector, and e is the random error 
vector. In K model: y = ^- Zu ^- The Q + K model 
can be written in a matrix form as: j = Qv + XyS + + e, 
where X, yS, Q and v are the same as those mentioned 
above, Z is the kinship matrix, and w is a vector of random 
genetic effects u~N (0, 2 K), where K is the kinship 
matrix. Manhattan plots were displayed using R software 
(v2.14.2). 

PGR amplification and sequencing 

Genomic DNA was extracted from leaves of 150 barley 
seedlings using CTAB method [29]. The primers used for 
DNA amplification were designed using the primer design 
tool of NCBI (http://www.ncbi.nlm.nih.gov/tools/primer- 
blast/). Primers for HvLDI were: forward, TTTTCGCA 
TGTCACCAAAAATGT; reverse, TCCGCTTCATTACCT 
TGGCG. The amplified DNA fragment was about 900 bp. 

The PCR reactions were completed as follows: 2.5 \A of 
10 X TransTaq HiPi buffer I (200 mM Tris-HCl (pH 8.4), 
100 mM (NH4)2S04, 20 mM MgCl2, 200 mM KCl, 2 \A of 
2.5 mM dNTPs, 1 \A of 10 \xM forward primers, 1 (il of 
10 (iM reverse primer, 0.5 \A of 5 units (iP^ of TransTaq 
polymerase High Fidelity (Beijing TransGen Biotech Co., 
Ltd.), and 1 \A of 50 ng of genomic DNA. The PCR ampli- 
fication program started at 5 min at 95°C, followed by 
32 cycles of 30 s at 95°C, 30 s at 60°C, and 1 min at 72°C, 
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Figure 1 The distribution frequency of LDI activity of the examined Tibetan wild barley accessions in 201 1 and 2012. The X axis 
represents LDI activity, the Y axis represents frequency of LDI activity. 



and then 10 min at 72°C for a final extension. The PGR 
products were analyzed by 1% agarose gel electrophor- 
esis. They were sequenced using 3730 x 1 DNA Analyzer 
(Applied Bio systems Inc., USA). 

Polymorphism of HvLDI and haplotype analysis 

Alignment and polymorphism detection of HvLDI se- 
quences were analyzed by MEGA v5.2 [30]. The haplotype 
analysis was conducted with Dnasp v5.0 [31]. 

Availability of supporting data 

The sequences of the 6 haplotypes were deposited in 
GenBank, and the accession numbers of these haplotypes, 
H1-H6 are KJ710426, KJ710427, KJ710427, KJ710429, 
KJ710430 and KJ710431, respectively. 

Results 

LDI activities of 162 Tibetan wild accessions 

A wide genetic variation of the LDI activity was observed 
among 162 Tibetan barley accessions for two years 
(Figure 1), ranged from -0.034 mU/(ig to 0.295 mU/(ig 
for 2011 and from 0.003 mU/(ig to 0.287 mU/(ig for 
2012, with a mean of 0.176 mU/(ig and standard differ- 
ence (SD) of 0.057 for 2011 samples, and a mean of 
0.167 mU/(ig and SD of 0.055 for 2012. 

Table 1 Variance component and heritabiiity (H^) 
estimates for LDI activity of 201 1 and 2012 





Vg 


Ve 




Year 201 1 


3.20E-03 


2.88E-04 


0.928 


Year 2012 


2.48E-03 


4.93 E-04 


0.854 



The results of variance analysis are presented in Table 1. 
Estimates of the variance component for genotype and 
error were 0.00320 and 0.000248 for 2011, 0.00288 and 
0.000493 for 2012. The heritabilities for 2011 and 2012 
were 0.928 and 0.854, respectively. 



Population structure 

The largest value of statistic index Ak was used as an indi- 
cator of the most probable number of subpopulations for 
all accessions (Figure 2). In the present study, the 162 
accessions could be classified into five subpopulations, 
with 19, 17, 21, 50 and 55 accessions for the each individ- 
ual subpopulation (Figure 3). The result was consistent 
with the data from the cluster analysis (Additional file 2: 
Figure SI). The population structure of 162 Tibetan wild 
barley accessions was listed in Additional file 3: Table S2. 



Note: Variance components are shown for Genotype (Vg), Error (Ve). 




Figure 2 Estimation of the most probable number of clusters 
(k), based on 10 independent runs and k ranging from 1 to 11. 
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Figure 3 Population structure of 162 Tibetan barley accessions. Population structure of 162 barley accessions was divided based on genetic 
diversity detected by 835 DArT markers with k = 5. Five subpopulations were represented by different colors. 



Genome-wide association study 

The LDI activities of 162 accessions and 771 DArT markers 
were used to perform genome-wide association study 
(GWAS). Two markers, bpb-8347 and bpb-0068 showed 
significance under Q, K and Q + K methods (Additional 
file 4: Table S3). When Q method was used, the -logio(p) 



values for bpb-8347 and bpb-0068 were 5.48 and 3.75 in 
2011 and 7.76 and 4.26 in 2012. If only kinship was con- 
sidered, the -logio(p) values for bpb-8347 and bpb-0068 
were 2.40 and 2.25 in 2011, and 3.10 and 2.15 in 2012. 
When a more restrict method, Q + K method was applied, 
the -logio(p) values were 2.75 and 2.30 in 2011, and 3.74 
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Figure 4 GWAS of LDI activity of 162 Tibetan barley accessions. GWAS was analyzed by three methods: a, d and g: Q method for 201 1, 2012 
and mean values of 201 1 and 2012; b, e, h: K method for 201 1, 2012 and mean values; b, d and f Q + K method for 201 1, 2012 and mean values. 
Significant associations were identified using criterion of-logio(P) >2 (P<0.01). 
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and 2.29 in 2012, all of which reached the threshold of sig- 
nificance in association analysis: -loglO(P) = 2 (Figure 4). 
For the mean LDI activity of the two years, the -logio(p) 
values for bpb-8347 and bpb-0068 were 8.07 and 4.88, 
respectively, if considered from only one population 
structure, and 3.06 and 2.66 if only kinship being con- 
sidered, and 3.70 and 2.71, if both population structure 
and kinship being considered. 

Bpb-8347 and bpb-0068, located at 73.44 cM on 
chromosome 6H and 66.50 cM on 3H, showed significant 
associations with LDI activity, explaining 6.8% and 5.3% of 
phenotypic variation in 2011 and 10.4% and 5.5% in 2012 
under Q + k method. Similarly, the two DArT markers 
could explain 10.0% and 6.6% of phenotypic variation for 
the mean values of 2011 and 2012. 



Association between HvLDI and LDI activity and 
haplotype analysis 

The variation of HvLDI and its promoter region were 
further investigated because several evidences supported 
the existence of a major gene responsible for LDI activity 
[7,10]. Alignment of this region revealed 35 SNPs in the 
150 accessions. Thirty-one SNPs (green-labeled in Figure 5) 
were significantly associated with LDI activity, explaining 
13.4% of the phenotypic variation and -logio(p) value of 
each SNP reaching a highly significant level (Figure 5). 

According to polymorphism of 35 SNPs, the 150 barley 
accessions can be divided into 6 haplotypes: H1-H6. 
H1-H4 showed high similarity of HvLDI sequence, while 
H5 and H6 had the similar sequences (Figure 5). Be- 
cause the 31 significant SNPs in HvLDI were complete 
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disequilibrium (r^ = 1) with each other, they could be 
regarded as one variant (Additional file 5: Figure S2). 
Hence, all accessions may be classified into two haplo- 
types, i.e. HA and HB. Accordingly, the distribution of 
haplotypes within different subpopulations was investi- 
gated. Ql, Q2 and Q5 were mainly constructed by haplo- 
type 1, and Q3 by haplotype 5, whereas Q4 subpopulation 
showed specificity to other subpopulations for its balanced 
composition of HA and HB (HA: HB = 18:29). Interestingly, 
HA showed significantly lower LDI activity than HB in Q4 
subpopulation (Figure 6b,c),while there was no difference 
between them in other subpopulations. 

Discussion 

In this study, we investigated the genetic variation of 
LDI activity and performed GWAS analysis of HvLDI, in 
order to identify the loci controlling LDI activity in Tibetan 
wild barley. A wide variation in LDI activity was found in 
the 162 Tibetan wild barley accessions, and some of which 
displayed a minus value. It should be noted that some 
dextrins released by the de-branching enzymes during ger- 
mination could also be detected by the Limit- Dextrizyme 
method, and these dextrins may account for the overesti- 
mation of LD activity [11,32,33]. Thus, the detected LD 
activities in extraction solutions of some accessions were 
higher than those of added LD, resulting in minus values 
of LDI activity, which was calculated as the reduced 
LD activity. However, the minus values of LDI in this 
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Figure 6 The distribution and LDI activity of different haplotypes. 

(a) The distribution of liaplotypes among different subpopulations. 

(b) Tlie LDI activity of HA and HB in total population and Q4 in 201 1 . 

(c) The LDI activity of HA and HB in total population and Q4 in 201 2 
HA includes HI, H2, H3 and H4; HB encompasses H5 and H6. Each bar 
is the mean of LDI activities of accessions belonging to the 
corresponding haplotypes. Error bars are SE values. 



experiment only occurred in very few accessions, and 
were very small (close to zero). Therefore the accessions 
with minus values of LDI are quite low in LDI activity. 
Actually, the impact of these dextrins produced during 
short-time extraction on the final LDI values could be 
negligible, as reported by [32]. 

Two DArT markers, i.e. bpb-8347 and bpb-0068 were 
found to be associated with LDI activities in Tibetan wild 
barley. Bpb-8347, a marker with highest GWA peak, is 
located near the locus of HvLDI gene obtained from IPK 
database, where the whole barley genome sequences lay 
[34]. This result confirms that HvLDI is a major gene con- 
trolling LDI activity in barley grain. Bpb-0068 might be a 
potential locus contributing to LDI activity, and it has not 
been reported up to date. Hence, it may be assumed that 
some minor effect genes affecting LDI activity may exist 
in Tibetan wild barley, and need to be confirmed by QTL 
analysis in segregating population in future work. 

A total of 31 SNPs of HvLDI were significantly associated 
with LDI activity. It is also found that 11 of these SNPs 
are significantly associated with the malt extract in 56 cul- 
tivated genotypes in a previous study [19], in which associ- 
ation analysis between HvLDI gene and LDI activity was 
not conducted. These results indicate that the LDI activity 
has a close relationship with malt extract, because the 
two traits are affected by the same 11 SNPs. Indeed, it 
was previously reported that the HvLDI antisense lines 
showed unpredicted pleiotropic effects on numerous 
enzyme activities, including a- and p-amylases, starch 
synthases, as well as starch granule types, all of which 
affect malt extract [35]. In addition, the LD activity- 
associated SNP in the previous study [19] was not detected 
in the Tibetan wild population, suggesting that LD activity 
in malt may be independent on polymorphism of HvLDI 
gene in Tibetan population. 

Interestingly, HA and HB showed a significant difference 
of the LDI activity only in Q4 subpopulation (Figure 6b,c). 
Q4 contains all kinds of haplotypes, while every other sub- 
population mainly consists of one haplotype (Figure 6a). 
In other word, Q4 has a wider genetic diversity than every 
other subpopulation. It seems that these associated SNPs 
are false positive due to highly structured LD (Additional 
file 5: Figure S2) in this chromosomal region, and similar 
phenomenon were also observed in association ana- 
lysis of frost tolerance [36,37]. It may be hypothesized 
that the variation of LDI activity in Q4 may be caused 
by other variations, which share linkage disequilibrium 
with the detected SNPs, in or around this region [34,38]. 
These variations should only belong to the accessions of 
HA in Q4 subpopulation, and may exist in the upstream 
or promoter region of HvLDI, or the regulated factors, 
which co-segregated with the detected SNPs. However, 
the assumption needs to be confirmed in the further 
work. 
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Four SNPs (13, 24, 72, and 352) within a coding region 
generate the alteration of the amino acids sequence. But 
these two kinds of haplotypes showed no significant 
difference in LDI activity among all accessions. It can be 
concluded that the alteration of LDI activity may be not 
due to the change of amino acid composition. 

Conclusions 

This study confirmed the locus of LDI gene and detected 
a DArT marker, bpb-0068, associated with LDI activity. 
The results provide useful information for identifying 
minor effect genes affecting LDI activity, and also prove 
that Tibetan wild barley is an elite germplasm, which 
may act as an abundant gene pool for barley breeding. 
Furthermore, the detected SNPs in this study should be 
helpful for better understanding the genetic control of 
LDI activity in barley grain. 

Additional files 



Competing interests 

The authors declare that they have r^o competir^g ir^terests. 
Authors' contributions 

YH, GZ and CL designed the experiment, YH and SC determined the LDI 
activity, GWAS analysis and drafted the manuscript. LY and YH participated in 
the sequence alignment. DW and FD performed the data analysis. All 
authors read and approved the final manuscript. 

Acknowledgments 

This research was supported by National Natural Science Foundation of 
China (Nos. 31330055, 31271656, 31129005), Zhejiang Provincial Natural 
Science Foundation of China (No. Z31 10054) and China Agriculture Research 
System. We thank Mr Xianyin Zhang and Ms. Mei Li, the technicians of 
985-lnstitute of Agrobiology for providing convenience in using the 
experimental equipments. 

Author details 

^Agronomy Department, Key Laboratory of Crop Germplasm Resource of 
Zhejiang Province, Zhejiang University, Hangzhou 310058, China. 
^Department of Agriculture and Food, Western Australia, WA 6983, Australia. 

Received: 27 January 2014 Accepted: 28 April 2014 
Published: 1 May 2014 

References 

1. Fincher GB: Molecular and cellular biology associated with endosperm 
mobilisation in germinating cereal grains. Anna Rev Plant Physiol Plant Mol 
6/0/1989, 40:305-346. 



2. Vasanthan T, Hoover R: Barley starch: production, properties, modification 

and uses. In Starch: Chemistry and Technology. 3rd edition. Edited by 
BeMiller J, Whistler R. New York: Academic press; 2009:601-628. 

3. Manners DJ: Starch degradation during malting mashing. Brewers Digest 
1974, 49:56-62. 

4. Manners DJ, Marshall JJ, Yellowlees D: The specificity of cereal limit 
dextrinase. Biochem J 1 970, 1 1 6:539-541 . 

5. Panchal CJ, Stewart GG: Utilisation of wort carbohydrates. Brewers Digest 
1979, 54:36-43. 

6. Sissons MJ, Lance RCM, Wallace W: Bound and free forms of barley limit 
dextrinase. Cereal Chem 1994, 71:520-521. 

7. Macri LJ, MacGregor AW, Shroeder SW, Bazin SL: Detection of a limit 
dextrinase inhibitor in barley. J Cereal Sci 1993, 18:103-106. 

8. MacGregor AW, Donald LJ, MacGregor EA, Duckworth HW: Stoichiometry 
of the complex formed by barley limit dextrinase with its endogenous 
inhibitor. Determination by electrospray time-of-flight mass spectrometry. 
J Cereal Sci 2003, 37:357-362. 

9. Sissons M, Taylor M, Proudlove M: Barley malt limit dextrinase: its 
extraction, heat stability and activity during malting and mashing. j/\m 
Soc Brew Chem 1 995, 53:1 04-1 1 0. 

10. Karakousis A, Barr AR, Kretschmer JM, Manning S, Logue SJ, Roumeliotis S, 
Collins HM, Chalmers KJ, Li CD, Lance RCM, Langridge P: Mapping and QTL 
analysis of the barley population Galleon x Haruna Nijo. Aust J Agr Res 
2003, 54(12):1 131-1 135. 

1 1 . Stahl Y, Alexander RD, Coates S, Bryce JH, Jenkinson HR, Morris PC: The 
barley limit dextrinase inhibitor: Gene expression, protein location and 
interaction with 14-3-3 protein. Plant Sci 2007, 172:452-461. 

12. Wu DZ, Qiu L, Xu LL, Ye LZ, Chen MX, Sun DF, Chen ZH, Zhang HT Jin XL, 
Dai F, Zhang GP: Genetic Variation of HvCBF Genes and Their Association 
with Salinity Tolerance in Tibetan Annual Wild Barley. PLoS One 201 1, 
6(7):e22938. 

13. Dai F, Nevo E, Wu D, Comadran J, Zhou M, Qiu L, Chen ZH: Avigdor 
Belles, Chen GX, Zhang GP: Tibet is one of the centers of 
domestication of cultivated barley. Proc Natl Acad Sci U S A20]2, 
109(42):1 6969-1 6973. 

14. Nguyen Viet L, Oene D, Marcos M, Benjamin K, Andreas G, Richard GF, 
Visser C, Gerard Van Der L: Association mapping of salt tolerance in 
barley (Hordeum vulgare L.). Theor AppI Genet 2013. doi:10.1007/ 
sOOl 22-01 3-21 39-0. 

15. Long NV, Dolstra 0, Malosetti M, Kilian B, Graner A, Visser RG, Van Der 
Linden CG: Association mapping of salt tolerance in barley (Hordeum 
vulgare L). Theor AppI Genet 2013. doi:l 0.1 007/sOOl 22-01 3-21 39-0. 

16. Neumann K, Kobiljski B, Dencic S, Varshney RK, Borner A: Genome-wide 
association mapping: a case study in bread wheat (Triticum aestivum L). 
Mol Breeding 2011, 27:37-58. 

17. Pasam RK, Sharma R, Malosetti M, Van Eeuwijk FA, Haseneyer G, Kilian B, 
Graner A: Genome-wide association studies for agronomical traits in a 
world wide spring barley collection. BMC Plant Biol 2012, 12:16. 

18. Cai SG, Yu G, Chen XH, Huang YC, Jiang XG, Zhang GP, Jin XL: Grain 
protein content variation and its association analysis in barley. BMC Plant 
Biol 2013, 13:35. 

19. Jin XL, Cai SG, Ye LZ, Chen ZH, Zhou MX, Zhang GP: Association of HvLDI 
with limit dextrinase activity and malt quality in barley. Biotechnol Lett 

2013, 35:639-645. 

20. Kristensen M, Planchot V, Abe J, Svensson B: Large-scale purification and 
characterization of barley limit dextrinase, a member of the alpha 
amylase structural family. Cereal Chem 1998, 75:473-479. 

21. MacGregor AW, Macri LJ, Schroeder SW, Bazin SL: Purification and 
characterization of limit dextrinase inhibitors from barley. J Cereal Sci 
1994 20:33-41. 

22. Pritchard JK, Stephens M, Donnelly P: Inference of population structure 
using multilocus genotype data. Genetics 2000, 155:945-959. 

23. Wenzl P, Li H, Carling J, Zhou M, Raman H, Paul E, Heamden P, Maier C, Xia 
L, Caig V, Jaroslava 0, Mehmet C, David P, Junping W, Rosy R, Smith KP, 
Muehlbauer GJ, Chalmers KJ, Andris K, Eric H, Andrzej K: A high-density 
consensus map of barley linking DArT markers to SSR, RFLP and STS loci 
and agricultural traits. BMC Genomics 2006, 7(1):206. 

24 Wenzl P, Carling J, Kudrna D, Jaccoud D, Huttner E, Kleinhofs A, Kilian A: 
Diversity Arrays Technology (DArT) for whole-genome profiling of barley. 
Proc Natl Acad Sci USA 2004 1 5(26):991 5-9920. doi:l 0.1 073/ 
pnas.0401076101. 



Additional file 1: Table SI. Minor allele frequency (MAF) and 
polymorphic information content (PIC) of 835 DArT markers used in this study. 

Additional file 2: Figure SI. Phylogenetic tree (neighbor-joining) of 
162 barley accessions based on 835 DArT markers. 

Additional file 3: Table S2. The value of population structure of 162 
accessions. Each accession belongs to the population with the highest 
value calculated by STRUCTURE software. 

Additional file 4: Table S3. GWAS results of bpb-8347 and bpb-0068 
under Q, K, Q + K methods. 

Additional file 5: Figure S2. Decay of linkage disequilibrium (LD) of 
HvLDI gene in Tibetan wild barley. Each point in the LD matrix represents 
a comparison between a pair of polymorphic sites. Different colors 
represent different levels of LD. 



Huang et al. BMC Plant Biology 2014, 14:1 1 7 
http://www.bionnedcentral.conn/1471-2229/14/117 



Page 8 of 8 



25. Alsop BP, Farre A, Wenzl P, Wang JM, Zhou MX, Romagosa I, Kilian A, 
Steffenson BJ: Development of wild barley-derived DArT markers and 
their integration into a barley consensus map. Mol Breeding 201 1, 
27:77-92. doi:1 0.1 007/s1 1 032-01 0-941 5-3. 

26. Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of 
individuals using the software STRUCTURE: a simulation study. Mol Ecol 
2005, 14:2611-2620. 

27. Hardy OJ, Vekemans X: Spagedi: a versatile computer program to analyze 
spatial genetic structure at the individual or population levels. Mol Ecol 
Notes 2002, 2:618-620. 

28. Bradbury PJ, Zhang ZW, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES: 
TASSEL: software for association mapping of complex traits in diverse 
samples. Bioinformotics 2007, 23(19):2633-2635. 

29. Huang JC, Ge XJ, Sun M: Modified CTAB protocol using a silica matrix for 
isolation of plant genomic DNA. Biotechniques 2000, 28:432-434. 

30. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGAS: 
Molecular evolutionary genetics analysis using maximum likelihood, 
evolutionary distance, and maximum parsimony methods. Mol Biol Evol 
2011,28:2731-2739. 

31. Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of 
DNA polymorphism data. Bioinform AppI Note 2009, 25:1451-1452. 

32. MacGregor AW, Bazin SL, Schroedor SW: Effect of starch hydrolysis 
products on the determination of limit dextrinase and limit dextrinase 
inhibitors in barley and malt. J Cereal Sci 2002, 35:17-28. 

33. McDougall GJ, Ross HA, Swanston JS, Davies HV: Limit dextrinase from 
germinating barley has endotransglycosylase activity, which explains its 
activation by maltodextrins. Plonto 2004, 218:542-551. 

34. Mayer KF, Waugh R, Brown JW, Schulman A, Langridge P, Platzer M, Fincher 
GB, Muehlbauer GJ, Sato K, Close TJ, Wise RP, Stein N, International Barley 
Genome Sequencing Consortium: A physical, genetic and functional 
sequence assembly of the barley genome. Nature 2012, 491:71 1-717. 

35. Stahl Y, Coates S, Bryce JH, Morris PC, Morris PC: Antisense downregulation 
of the barley limit dextrinase inhibitor modulates starch granule size 
distribution, starch composition and amylopectin structure. Plant J 2004, 
39:599-611. 

36. Fricano A, Rizza F, Faccioli P, Pagani D, Pavan P, Stella A, Rossini L, Piffanelli 
P, Cattivelli L: Genetic variants of HvCbf14 are statistically associated with 
frost tolerance in a European germplasm collection of Hordeum vulgare. 

Theor AppI Genet 2009, 1 1 9:1 335-1 348. 

37. Visioni A, Tondelli A, Francia E, Pswarayi A, Malosetti M, Russell J, Thomas W, 
Waugh R, Pecchioni N, Romagosa I, Jordi C: Genome-wide association 
mapping of frost tolerance in barley (Hordeum vulgare L.). BMC Genomics 
2013, 14:424. 

38. Myles S, Davison D, Barrett J, Stoneking M, Timpson N: World wide 
population differentiation at disease associated SNPs. BMC Genonnics 
2008, 1:22. 



doi:1 0.1 1 86/1 471 -2229-1 4-1 1 7 

Cite this article as: Huang et al.: Genetic architecture of limit dextrinase 
inhibitor (LDI) activity in Tibetan wild barley. BMC Plant Biology 
2014 14117. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at /^\ Ri^nHod Tpntral 

www.biomedcentraLcom/submit momea L.enTrai 



